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Abstract: Based on the magnitude-redshift diagram for the sample of supernovae la analysed by Perl- 
mutter et al. (1999), Davis & Lineweaver rule out the special relativistic interpretation of cosmological 
redshifts at a confidence level of 23 (t. Here, we critically reassess this result. Special relativity is known 
to describe the dynamics of an empty universe, by means of the Milne kinematic model. Applying 
only special-relativistic concepts, we derive the angular diameter distance and the luminosity distance 
in the Milne model. In particular, in this model we do not use the underlying metric in its Robertson- 
Walker form, so our exposition is useful for readers without any knowledge of general relativity. We 
do however, explicitly use the special-relativistic Doppler formula for redshift. We apply the derived 
luminosity distance to the magnitude-redshift diagram for supernovae la of Perlmutter et al. (1999) 
and show that special relativity fits the data much better than that claimed by Davis & Lineweaver. 
Specifically, using these data alone, the Milne model is ruled out only at a 2a level. Although not 
a viable cosmological model, in the context of current research on supernovae la it remains a useful 
reference model when comparing predictions of various cosmological models. 
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1 Introduction 

In a recent paper. IPavis fc Linewea^^ (l2004fl attempt¬ 
ed to clarify several common misconceptions about the 
expansion of the universe. In particular, they convinc¬ 
ingly pointed out that uniform expansion of an infinite 
universe implies that very distant galaxies recede from 
us with superluminal (faster-than-light) recession ve¬ 
locities. Moreover, we can observe such galaxies. This 
does not violate special relativity (SR), because their 
velocities are not measured in any observer’s inertial 
frame. They are measured in the so-called reference 
frame of Fundamental Observers, for which the uni¬ 
verse looks homogeneous and isotropic. 

Unfortunately, Davis & Lineweaver not only clar¬ 
ified some misconceptions; they also created a new 
misunderstanding. They claimed to “observationally 
rule out the SR Doppler interpretation of cosmological 
redshifts at a confidence level of 23(T.” (The special- 
relativistic interpretation of redshift is kinematic, i.e., 
the Doppler effect.) As we will explain later, Davis 
& Lineweaver did not apply SR properly. Specifically, 
they did not consistently use the definition of the lu¬ 
minosity distance, Dl- SR is known to describe the 
dyna mics of an empty universe iPeacoc kllQQfltlLongairl 
I 2 OO 3 II . Their error led Davis & Lineweaver to an ex¬ 
pression for the luminosity distance as a function of 
redshift that was entirely different from that for an 
empty universe. 

In the framework of general relativity, the calcu¬ 
lation of Dl for an empty universe is straightforward. 
However, here we will present also an alternative ap¬ 
proach, based entirely on SR. Namely, we will de- 
ri ve Dl applyi ng the kinematic cosmological model 
of lMilnd (Il933fl . This approach is useful for readers 


without any knowledge of general relativity. It eluci¬ 
dates the meaning of time and kinematics in cosmol¬ 
ogy. The Milne model offers also an interesting insight 
into all Friedman-Robertson-Walker (FRW) cosmolog¬ 
ical models. 

The paper is organized as follows. We begin in 
Sec. H by deriving the angular diameter distance for 
an empty universe using the FRW framework. Next, 
we present the corresponding derivation in the Milne 
model. In Sec. Owe derive the luminosity distance in 
the Milne model and compare it to the angular diam¬ 
eter distance. Finally, in Sec. Owe present the result¬ 
ing magnitude-redshift diagram for supernovae la. We 
summarize in Sec.O 

2 Angular diameter distance 

For comparison, let us first recall the derivation of the 
angular diameter distance for an empty universe in the 
FRW framework. The metric of a homogeneous and 
isotropic universe is given by the Robertson-Walker 
line element: 

c?ds^ = c?dt^ — a^{t)[dx^ -\- {x/Ro)dD\. (1) 
Here, 

dO. = d6^ sin^ 6d(f>^, (2) 

and Ro^ is the (present) curvature of the universe. 
The function S{x) equals sin( 2 :), x, and sinh(a:) for 
a closed, flat, and open universe, respectively. The 
function a{t) is called a scale factor and relates the 
physical, or proper, coordinates of a galaxy, r, to its 
fixed or comoving coordinates, x: r = ax. This func¬ 
tion accounts for the expansion of the universe; its 
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detailed time dependence is determined by the Fried¬ 
man eqnations. We normalize a so that at the present 
time, a{to) = 1. For an open universe, the Friedman 
equations yield 

cH-^ 

° (l-no)l/2’ 

where Ho is the present value of the Hubble constant 
and Ho is the present value of the total (energy) den¬ 
sity in the universe, in units of the so-called critical 
(energy) density. 

Let us work out the angular size of an object of 
proper length Ay, perpendicular to the radial coordi¬ 
nate at redshift 2 . The relevant spatial component of 
the metric 0-0 is the term in dO. The proper length 
Ay of an object at redshift z, corresponding to scale 
factor a(te), for an open universe is 

Ay = a{te)Ro sinh(a:/i?o)A0 = DaA6, (4) 


where we have introduced the angular diameter dis¬ 
tance Da ~ a{ts)RoSmh{x/Ro)- Here, te is the time 
of emission of photons. Since a{te)~^ = 1 -|- 2 , we have 


Da = Rosmh{x/Ro)/{l + z). (5) 

For an empty universe. Ho = 0, hence Ro = cHo^. 
Moreover, then a{t) — Hot and the equations of null 
radial geodesics are easy to integrate. The result is a; = 
cHo^ ln(l-|- 2 ). Substituting this into equation 0 and 
using the definition of the hyperbolic sine, we obtain 


Da{z) = cHo^ 


z{l + zl2) 
(l + z)^ 


( 6 ) 


This is the angular diameter distance for an empty 
universe, derived in the FRW framework. 

In the Milne model, the cosmic arena of physi¬ 
cal events is the pre-existing Minkowski spacetime. In 
the origin of the coordinate system, O, at time t = 0 
an ‘explosion’ takes place, sending radially fundamen¬ 
tal observers with constant velocities in the range of 
speeds (0, c). The fundamental observer with velocity 
V, Fo, carries a rigid rod of length Ay, oriented per¬ 
pendicularly to the line of sight of the observer at O. 
At time to this rod emits photons. At the photons’ 
arrival time at O, to, the rod subtends at O an angle 


Acj) = Ay/ro, (7) 


where re is the distance from O to at the time of 
emission of the photons, to- We have to = to + U, 
where tt is the travel time of the photons. Since 



Te 

te = — 

(8) 


V 

and 

tt = —, 

(9) 


c 

we obtain 

r -ct ^ 

(10) 


where [3 = v/c. The special-relativistic formula for the 
Doppler effect is 


where 2 is the photons’ redshift; hence. 


(1 + 2 )^ 
^ {1+zr + i- 

Using equation in Gnj yields 


re = Cto 


z{l + z/2) 

(l + z)^ 


( 12 ) 


(13) 


Since in the Milne model, for any time t and for any 
fundamental observer F^, r = vt, the observer at O 
observes the Hubble flow: v = Hr, where the Hubble 
constant is H = t~^. Hence, to = HF^- The angular 
diameter distance is defined via the equation Ay = 
DAA(j>. Therefore, using equation 0 we obtain finally 


Da{z) = re = cHo ^ 


2(1-b 2/2) 
(l + z)^ 


(14) 


In the above derivation we have applied special rel¬ 
ativity a number of times. First, writing Eq. 0 we 
have followed its central assumption, that the velocity 
of light is always c, regardless the relative motion of the 
emitter and the observer. Secondly, we have applied 
the special-relativistic, kinematic interpretation of red¬ 
shift - i.e., the Doppler effect - and used formula o 
for it. Finally, writing Eq. 0 we have assumed that 
the geometry of space is Euclidean. 

Eq. II41 exactly coincides with formula 0 for Da 
for an empty universe, the latter derived from the met¬ 
ric in its FRW form. Our second derivation, which 
was purely special-relativistic, employed other (i.e., 
conventional Minkowskian) dehnitions of distance and 
time. Still, we arrived at the same formula for Da- 
This is so because this formula relates direct observ¬ 
ables: proper (i.e., rest-frame) size of an object to its 
angular size and redshift. Regardless of what are the 
dehnitions of coordinates in a given coordinate system, 
their consistent application should lead to the same re¬ 
sult in terms of observables! 


3 Luminosity distance 

The relation between the luminosity distance, Dl , and 
the angular diameter distance is 

DL = {l + zfDA. (15) 

However, this is true in general relativity; let’s check 
whether this holds also in the Milne model. Let’s place 
a source of radiation at the origin of the coordinate 
system, O. We assume again that at time t = 0 at O an 
‘explosion’ takes place, sending radially fundamental 
observers with constant velocities in the range of (0, c). 
At time te the source emits photons, which at time to 
reach a fundamental observer moving with velocity v, 
such that 

vto = c{to — te). (16) 

If the source emits continuously photons with constant 
bolometric luminosity L, the observer at To = vto re¬ 
ceives a flux of radiation with bolometric intensity 



1 -b2 = 


1-/3 


47rri(l + 2 )^ 


(17) 










www.publish.csiro.au/joumals/pasa 


3 


The factor (1 + z)^ in the denominator is due to the 
Doppler effect. Specifically, one factor 1 + z is due 
to the fact that the wavelength, and so the energy, of 
the observed photons is redshifted. The second factor 
1 + 2 is due to the fact that photons, emitted in the 
time interval Ate, arrive to the observer in the time 
interval to = {1 + z)Ate. 

The luminosity distance is defined by the equation 
/ = Z//(47rD|,), hence 

Dl = {1 + z)ro = {1 + z)cl3to. (18) 

Time to is the time of observation indicated by the 
clock at the source, but the fundamental observer is 
moving with respect to O, so according to SR, his 
clock delays compared to that at O. At the moment 
of observation, his clock shows time To = to/^, where 
7 = (1 — Therefore, we have 

Dl = {1 + z)cTo(3'y. (19) 

Relative to the observer, however, this is the source 
that is moving. Using his clock, the observer deduces 
that since the Big-Bang, the source has moved off to 
the distance r'o = vto- Because the source is mov¬ 
ing, the distance r'o is length-contracted relative to its 
rest-frame value, Vo'. r'o = ro/y. Hence, the observer 
will agree that To = yvTo = CTofiy, what yields equa¬ 
tion 11911 . 

Next, we have 

^ (1 - /3)i/2(i + /3)i/2 = 


relative to O and R denotes the position of O' relative 
to O. Hence, v' = Ho{r — R) = Hov' , where r' is the 
position of O" relative to O'. Thus, the observer O' 
observes an isotropic Hubble flow around him, so he 
is apparently at the center of expansion. The point 
is that this r esult holds al so for relativistic velocities 
llMilnelll93.3l : lR,indlerlll977T> . Strictly speaking, in the 
Hubble law the Hubble constant is an inverse of the 
local proper time of the observer at O'. This is why in 
Eq. iU we have identified Ho with 

It is rather surprising why Milne, who insisted so 
much on an observables-oriented approach to cosmol¬ 
ogy, did not derive himself an explicit formula for the 
luminosity distance as a function of redshift. We will 
see, however, that he was close to it. The Ap pendix 
to hi s classical paper on kinematic relativity dMii^ 
ll 9.3.311 bears the title “The apparent brightness of a re¬ 
ceding nebula”. Its final formula (16) describes “the 
total light received by A (the observer) on his own 
photographic plate” in a fairly complex, integral form. 
This formula involves implicitly the luminosity dis¬ 
tance, but in terms of the recession velocity of the neb¬ 
ula rather than its redshift. However, while the red- 
shift of a distant nebula is a direct observable, its ve¬ 
locity is not. Extracting the luminosity distance from 
the formula and using the SR relation between veloc¬ 
ity and redshift, Eq. m, we rederive the formula for 
the luminosity distance given by our Eq. (HIJ. 

4 Apparent magnitude — red¬ 
shift relation 


where in the last equality we have used the SR for¬ 
mula for redshift, Eq. GU- Combined with Eq. GU, 
Eqs. 119t - ll2()l l yield 

Dn = (1 + + = cHo^z{l + z/2). (21) 

(1 -I- z)-^ 

Comparing with Eq. 11411 we see that indeed Dl = 
(1 -1- z)^Da, in accordance with general relativity. 

To derive the angular diameter distance, in Sec.Gl 
we have used the observer’s rest-frame. To derive the 
luminosity distance, however, in the present section 
we have switched to the source’s rest-frame. We have 
done so for simplicity of the resulting calculations. In 
particular, only in the latter frame is the radiation of 
the source isotropic, and one can apply simple equa¬ 
tion m for the observed flux. 

Deriving the two distances, we have placed either 
the observer or the source at a special position: at 
the center of expansion. Are then our results gen¬ 
eral? Yes: although in the Milne model this center 
does indeed exist, every fundamental observer consid¬ 
ers himself to be at the center of expansion! This can 
be easily seen in the non-relativistic regime. According 
to the Galilean transformation of velocities, the veloc¬ 
ity of any observer O" relative to another observer O' 
is v' = V — Vrei, where v and Vrei are respectively the 
velocity of O" relative to O and the velocity of O' rel¬ 
ative to O. But by the construction of the model, the 
observer at O observes the Hubble flow, so v = Hof 
and Vrei = TIoR, where r denotes the position of O" 


The apparent bolometric magnitude, ms, of a stan¬ 
dard candle located at redshift 2 is related to its abso¬ 
lute bolometric magnitude, M, by the equation 

ms = 5 logio Dl + 25 + M. (22) 

Here, Dl is the luminosity distance, expressed in me¬ 
gaparsecs (Mpc). Supernovae la turned out to be very 
goo d standard candle s i n the universe (se e , for e xam¬ 
ple. iPerlmutteJ l2003^ . iPerlmutter et al.l lll99ltl cast 
the above equation to the form 


ms = 51ogiQ(HoDs)-I-A4, (23) 

where 

M = M - 5 logio Ho + 25 (24) 

is the magnitude zero-point, and the Hubble constant 
is expressed in km ■ s“^ • Mpc“^ (Eq. (l)-(2) of Perl- 
mutter et al. 1997). We prefer to rewrite Eq. 12311 to 
the form 


ms = 5 logio 



+ M, 


(25) 


where 


M = M + 5 logio c = M+ 25+ 5 logio 2-998; (26) 

the argument of logarithm in Eq. 1251 is explicitly di¬ 
mensionless. Introducing Eq. GB in 1251 yields 


ms = 5logio [-^(1 + '*/2)l + -M- 


(27) 
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This is the magnitude-redshift relation in the Milne 
model. It coincides exactly with the corresponding 
relation for an empty universe. In Fig.^it is shown as 
a dotted line. For the magnitude ze ro-point we adopt 
the va lue A4 = —3.32, calculate d in iPerlmutte r et alJ 
ill 99711 and used also in iPerlmntter et" .l lll999^ . 



Figure 1: The magnitude-redshift diagram for su¬ 
pernovae la. The thick solid line shows the pre¬ 
diction of the currently favored cosmological model, 
=0.28 and Ha = 0.72. The thin solid line shows 
the prediction of special relativity according to Davis 
& Lineweaver. The dotted line shows the correct pre¬ 
diction of special relativity. 


To derive the luminosity distance Dl = (1 -I- z)D, 
Davis & Lineweaver nsed the Hubble law: 


D = cH-^p. (28) 

Comparing with Eq. m we see that this leads to 
the correct expression for Dl provided we identify H 
with Ho‘‘\ i.e. the value of the Hubble constant at 
the source and at the time of observation of photons. 
However, they identihed H with i.e. the Hubble 
constant at the source but at the time of emission of 
photons: 

=c[Hi‘’^]~^ P = cPt,. (29) 

Combined with the SR expression for f3 as a function 
of redshift, Eq. d- the above equation yields 


j^{DL) ^ c (1 + g)^ - 1 
( 1 - 1 - 2)2 - 1 - 1 ' 


(30) 


(Eq. (10) of iDavis fc Lineweaveill2004ll . Next, Davis 
& Lineweaver applied the equality 


= {1 + z)Hi°\ (31) 


correct in case of an empty universe. Let’s check whe¬ 
ther this holds also in the Milne model. First, 11^“^ = 
and where t and r are the time mea¬ 

sured (by Fundamental Observers) respectively at the 
source and at the observing point. Hence, 




To -1 


— = 7 P 

Pe 


-1 


(32) 


where p = t^jto. From Eq. II til we have 
/3 = 1 - p, 
so 

7-' = (l-/3")'/" = b(2-p)ll/^ 

Used in Eq. iTTa . the above equation yields 

(s) /„ ..\l/2 


77, 


77, 


(o) 


2-p 

P 


In turn, Eq. lIMl . used in Eq. IlH . yields 

1/2 


l-b 2 = 


2-p 


(33) 

(34) 


(35) 


(36) 


By inspection of Eqs. igg-® we see that they in¬ 
deed imply Eq. m- Once again, we have found that 
the Milne model exactly describes the dynamics of an 
empty universe. 

Comparing equation with ilTSl we see that 
equation l|^ leads to an expression for the luminosity 
distance which underestimates the correct one by the 
factor of to/t^. We have to/te = {tolro) ■ (vo/te) = 
7(1 + z) = [(1 -|- z)^ + l]/2. Using Eq. 13H in Eq. 1301 
we obtain 


jj(DL) 




c {1 + zf-l 
77(0 {l + z)^ + l’ 


therefore 


d(dl) ^ = 


c (l + ^)^-l 

Ho {l + zy + P 


(37) 


(38) 


Here, Ho denotes the value of the Hubble constant 
at the observing point at the time of observation and 
for simplicity we have skipped the superscript “Oi’y 
Combined with Eq. ll^ . Eq. ll^ yields the function 
shown in Fig. Qas a thin solid line. Davis & Lineweaver 
did not provide an explicit formula for the magnitude- 
redshift relation in SR, but this line closely follows 
their plotted curve. 

The thick solid line in Fig.^shows the magnitude- 
redshift relation for the currently favored cosmological 
model: a flat universe with a nonzero cosmological con¬ 
stant, D„i = 0.28 and Da = 0.72.*^ Useful expressions 
for the lumin o sity d istance in this model are given in 
IChodorowskil (1200, *111 . We see that the thin solid line 
(prediction of SR according to Davis & Lineweaver) 
is very distant fro m the thick solid line. Using the 
supernovae data of IPerlmutter et all (Il99fltl . Davis & 
Lineweaver verihed that the model given by the thin 
line “is ruled out at more than 23 (j” compared with 
the currently favored model. 

On the other hand, the dotted line (correct pre¬ 
diction of SR) follows the t hick solid line much mor e 
closely. Given the data by IPerlmutter et a.ll ill 999^ , 

^ These are recent estimates of the cosmological parame¬ 
ters flm and Ua from the magnitude-redshift relation of 
observed high-redshift supernovae, combined with other 
observ ational constraints. See, for example, iTonrv et al.l 

|20 o3 ). 
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how close is the prediction of SR to that of the fa¬ 
vored model? The answer is provided by Perlmut- 
ter et al. themselves: their analysis yielded the con¬ 
straint O.Sflm — O.GflA = —0.2 ± 0.1. An empty uni¬ 
verse corresponds to Q.m = Ha = 0, hence for SR, 
O.Sflm, — O.GUa = 0. We see thns that within two 
standard deviations, the data was consistent with an 
empty universe! This is not to defend this model as 
a viable alternative to the currently favored model. 
From our mere existence we know the universe is not 
empty. A host of observational evidence consistently 
points towards the currently favored model. This is 
only to say that a few years ago, the evidence from su¬ 
pernovae data alone for the accelerated expansion of 
the universe was not so strong, and the assumption of 
its purely kinematic expansion at low redshifts could 
then serve as a reasonable starting approximation. 

Since the year 1999 the supernovae la data has 
impro ved. In particular, the analysis of iTonrv et alJ 
11200,^1 yielded Ha - 1.40^ = 0.35 ± 0.14. This con¬ 
stituted a modest improvement over the result of Perl- 
mutter et ah, implying the present acceleration of the 
universe’s expansion to be detected at a 2.5cr confi¬ 
dence.^ A significant improvement was achieved by 
discovering and observing supernovae la a t z > 1 with 
the H ubble Space Telescope. Fig. 8 of iRiess et a.lJ 
1 I 2 OO 4 II shows the resulting joint confidence intervals 
for (fim, flv) from SNe la. In this figure, the point 
Qm = Ha = 0 lies outs ide shown confidence contours 
of 99.7%. iRiess et al.l claim that “with the current 
sample, the 4a confidence intervals (i.e., > 99.99% 
confidence) are now fully contained within the region 
where Ha > 0.” Similar results of the analysis of the 
IRiess et alJ sample have been obtained independently 
w lWrightlll200^ . 

An empty universe is thus not a viable cosmologi¬ 
cal model, but remains a useful reference model when 
compari ng predictio ns of various cosmological models. 
Fig. 7 of iR.iess et al.l shows the magnitude-redshift dia¬ 
gram for SN la in a residual form: relative to an empty 
universe model. Being eternally coasting, this model 
has a vanishing deceleration parameter, so it naturally 
separates accelerating fro m deceleratin g models. Also, 
it is evident from Fig. 7 of lR,iess et alJ that the model 
fits the data much better than an Einstein-de Sitter 
universe {Qm = 1, Ha = 0). 


verse in the FRW cosmology. We have shown where 
Davis & Lineweaver failed to correctly derive the lu¬ 
minosity distance. Finally, we have presented the re¬ 
sulting magnitude-redshift diagram for supernovae la. 
While the prediction of special relativity according to 
Davis & Lineweaver is far away from that for the cur¬ 
rently favored cosmological model, the correct predic¬ 
tion of special relativity follows the favored model much 
more closely. Though not a viable alternative to the 
currently favored model, the Milne model has great 
pedagogical value, elucidating the kinematic aspect of 
the universe’s expansion. In the context of current re¬ 
search on supernovae la, it remains a useful reference 
model when comparing predictions of various cosmo¬ 
logical models. 
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5 Summary 

We have derived the angular diameter distance Da 
and the luminosity distance Dl in the Milne kinematic 
cosmological model, using only special-relativistic con¬ 
cepts. In the derivations, the central role was played 
by the special-relativistic Doppler formula for photons’ 
redshift. We have found that Dl = (1 z)^Da, in ac¬ 
cordance with general relativity. The derived formulae 
are identical to these corresponding to an empty uni- 

^From supernovae alone. Combined with the constraint 
of a flat universe, strongl^_su£ported by the CMB obser¬ 
vations, the data of |Tonr^_e^_aJJ yielded Qm = 0.28 ± 0.05 
and Da = 0.72 ± 0.05, implying a currently accelerating 
universe at much higher confidence. 


















